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We theoretically demonstrate a mechanism for bright soliton generation in spinor non-equilibrium 
Bose-Einstein condensates made of atoms or quasiparticles such as polaritons in semiconductor 
microcavities. We give analytical expressions for bright (half) solitons as minimizing functions 
of a generalized non-conservative Lagrangian elucidating the unique features of inter and intra¬ 
competition in non-equilibrium systems. The analytical results are supported by a detailed numerical 
analysis that further shows the rich soliton dynamics inferred by their instability and mutual cross¬ 
interactions. 
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INTRODUCTION 

Investigations on Bose-Einstein condensation (BEC) have been key in understanding quantum matter at very low 
temperature highlighting its distinctive coherent wave properties [HE]. BEC is a macroscopic quantum phenomenon 
often involving a large number of interacting coherent quantum particles [2]. In the last two decades the phase 
transition of many-body systems to a BE condensed state has become feasible within a great variety of different 
physical systems such as solid-state light-matter systems, various species of atoms and molecules and for gases of 
photons or classical waves GSfia. Among the most intriguing properties of all these condensed coherent matter states 
are regimes of frictionless flow (superfluidity) [7| below a critical velocity [H S3 or the response to motion via elementary 
excitations such as quantized vortex rings, quantum vortices and lattices thereof, dark and bright solitons [SHE]. 
Different kinds of BEC inherit properties depending on the nature of the entities constituting them, such as variable 
inter-particle interactions due to Feshbaclr resonances when applying certain external fields in atomic and polariton 
condensates [nnii] or local particle sources and sinks in non-equilibrium condensates of exciton-polaritons Bum 
or atoms [20] . In turn the out of equilibrium aspect of open quantum systems as well as the nature of interactions 
between the coherent particles is key for the possible pattern formations in those condensates and thus the properties 
of excitations within the condensate HII US- 

Attractive Bose-Einstein condensates of atoms were experimentally realized [21] as well as repulsive condensates of 
polaritons [5] . Later it was shown that self-interactions of condensed atoms or polaritons can be tuned continuously via 
utilizing Feshbach resonances mm providing one possible route towards implementing attractive non-equilibrium 
condensates. In addition quasiparticles in solid state systems can have an effectively negative mass proportional to 
the inverse of the second derivative of their dispersion, for which polaritons are an example, and so the mathematical 
form of the mean held theory corresponds to effectively attractive self-interactions providing an alternative approach 
to Feshbach resonances [221 123] . Polaritons are quasiparticles within semiconductor microcavities which satisfy Bose- 
Einstein statistics 0] and show macroscopic occupation of the lowest energy state with off-diagonal long range order 
[5] . Temporal coherence is limited to a few ps due to the continuous loss of particles and information [24] . Excitons 
are coupled electron-hole pairs of oppositely charged spin-half particles in a semiconductor bound by the Coulomb 
interactions m • Since states of excitons interact with light [25] in the strong coupling regime exciton-polaritons 
(polaritons) are formed [2(5] and as polaritons are 10 9 times lighter than rubidium atoms [5], condensation can be 
observed on the Kelvin range in CdTe/CdMgTe/GaAs micro-cavities [U (5( [27] and even at room temperatures in 
GaN or flexible polymer-filled micro-cavities [28l 129] . Polaritons are a prime example for an open quantum system 
and allow to test non-equilibrium properties and are analoguous in their mathematical structure to open atomic BEC 
systems when considering mean held properties of their coherent ground states [20 . 

The system we shall study in this paper is an effectively attractive spinor non-equilibrium condensate within a mean- 
held parameter regime and with potential realization in non-equilibrium systems of atoms, molecules or quasiparticles 
such as polaritons. A key requirement for the applicability of our theory is the continuous exchange of particles with an 
external bath, which allows richer dynamical behavior of the spinor held and the emergent excitations within. Effective 
mechanisms to generate excitations in general BEC within the Gross-Pitaevskii (GP) mean-held parameter regime 
Bill] are important tools for investigations on their properties both experimentally and theoretically [30] and carry 
the possibility for applications e.g. in optical computing nm EH- Vortices with quantized circulation in superhuids 
are a key example of a topologically stable excitation in a 2d BEC system guided by a nonlinear Schrodinger-type 


2 


equation. Recently, they have been predicted [52HM| and observed [35] in polariton condensates. Their counterparts 
in Id repulsive condensates are solitons, localized excitations moving without deformation over time |15j . They have 
been predicted in several contexts [HJ [361442] and observed so far mostly in planar cavities (22j [43| 04] . In addition 
state-of-the-art investigations describe driven dissipative bright solitons under non-resonant excitation [35] SB] for the 
repulsive BE condensate. In this work we provide an analytical mean field non-equilibrium theory for bright solitons 
and predict numerically the bright soliton patterns that emerge within an effectively attractive spinor condensate. 
Our numerical results are supported by analytical parameter equations and particularly explicit expressions for the 
stationary solutions are obtained by this method. 

The paper is organized as follows. First we introduce the state-of-the-art mean field model. Then we simplify the 
mean-field model to the Keeling/Berloff form of the governing non-conservative partial differential equation (PDE). 
Based on this we introduce the Lagrangian treatment of the PDE for a bright soliton solution. Subsequently we 
extend the considerations to the spin sensitive case. Finally we turn to numerical investigation of the bright solitons 
and particularly compare numerics with analytics. 


THE MODEL 

We now turn to the state equation for a two component spinor field ip = (ip+,ip-) T that describes the two coherent 
macroscopically occupied many body states ip± where we refer to ± as spin although ± can be thought of as any 
two coherent hyperhne states in the considered many body system. The system/scenario of interest is modeled by 
ID Ginzburg-Landau type equations coupled to a rate equation for the reservoir hr capturing the non-condensed 
particles within the bath feasible in open systems such as polaritons PG2 or atoms {2 ED], 

id t ipa = ( 1 ) 


with er £ {+, —} and 

= (1 - + ai(l ^±| 2 + n R) + « 2 |- 0 =fI 2 +U S jip± + ^ ('yriR - T) ip± ( 2 ) 

Here g is the amount of energy relaxation of the condensed phase, aq denotes the same spin interaction strength and 
«2 the cross spin interactions, m is the effective mass of the particles within the condensate phase. Furthermore we 
denote the pumping rate into the condensate 7 and the decay rate for condensed particles as F. Pumping is modeled 
by a spatial and time dependent function P which induces particles into the reservoir/bath of the system 15, 20] . 

d t n R = P - n R (T R + j{\ip + \ 2 + \ip„\ 2 )) (3) 

where [201 has been extrapolated to the spin sensitive case. The reservoir balance equation is fixed at a decay rate 
T^ and driven by the spatially dependent pump term P = Ap exp(— x 2 /a^) where ap is the width, and Ap is the 
amplitude modeling the scattering into the ground state. The reservoir scatters particles into the condensate at a 
rate 7 . The coupled system of equations describes the time evolution of the condensate and its reservoir and is an 
example for a non-resonant pumping scheme in polariton condensates Hang with analog in atoms [20]. For atoms 
one can understand the spinor field as the two hyperfine states, say A and B, of a BEC ip [2]. Now, let us turn to 
analytical considerations for the spin coherent case. 


ATTRACTIVE NON-EQUILIBRIUM CONDENSATE 

Let us first consider a simple case to introduce the analytical arguments presented in the following. In the case of 
spin coherence we can simplify the system of equations to [all [48] 

-d t ip = (1 — + a' 1 (\ip\ 2 +n R ) + Ujip+ - ('yn R ~T)ip. (4) 

In addition we note that the true form of the non-equilibrium dynamics has not yet been measured experimentally and 
both models [Ml 35] are often used to successfully describe experiments with polariton condensates. The actual form of 
the density saturation in Q however, becomes relevant for large amplitudes of the wave function, i.e. large occupation 
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numbers of the condensed phase. Both phenomenological mean field models are connected by Taylor expanding the 
pumping term stated in [48j to the first order as outlined in |ll 6 j . Now, starting from ([2]) we approximate/simplify the 
complex terms of the condensate wave equation Q by the reservoir and decay of particles as in (16] 1391 , 


l - ('yn M - T) ~ i(S + £a'i|i/>| 2 ), 


(5) 


to obtain a mathematical form easily accessible to apply the following analytics. For attractive condensates the 
coefficient of the real-valued nonlinear term a! x in Q is negative, while it is positive when modeling repulsively 
interacting condensates. Bright soliton are solutions for the GPE equation with negative nonlinearity and thus 
excitations in attractive BE condensates. In (49, 50 1 , [521 the bright soliton solutions for the complex cubic-quintic 
Ginzburg-Landau equation have been considered. By taking 0 into account eq. @ can be regarded as a special 
case of this cubic-quintic PDE, but without quintic terms and we can apply and later extend the analysis outlined in 
US [50) 152] . However, in contrast to their derivation we assume that the sign of the linear complex term is positive, 
while stability requires us to introduce a cutoff |34| and to set in turn e = — |e| which is the magnitude of the density 
dependent decay/saturation of condensed particles. Further we neglect the blueshift from the reservoir (which can be 
absorbed into the phase) and consider no external potential, i.e. U = 0 in ©• Rescaling r = yj2m / h ■ x, z = Kt and 
if) —> i— E yields the governing PDE for the non-equilibrium system, 


id z E= -AE± \E\ 2 E + i5E + ie\E\ 2 E ( 6 ) 

with the d dimensional Laplacian A E = r 1 _d i d r (r d ~ 1 d r E). In the following we assume d = 1 as we are interested 
in the description of bright solitons, but write down the Lagrangian procedure for general dimensions. For negative 
self-interactions 0 with 1/2 in front of the Laplacian and with <5 = 0 and e = 0 the equation possesses an exact 
bright soliton solution 0154] 


E(z, t) = A - sech(A(,z — vt)) exp (i{nz — wf)), (7) 

where A denotes the amplitude of the excitation, k the soliton wave number ui = ( k 2 — A 2 )/2 is the soliton dispersion 
relation and v = du/dn = k the soliton velocity . 541 . Later we extend 0 to the variational ansatz we will use to 
describe non-equilibrium bright solitons. The Lagrangian of our system Q denoted L, can be separated additively into 
its conservative part Lq (real-valued energy terms) and its dissipative part, inducing the non-equilibrium properties 
of the condensate wave function due to pumping and loss of the particles, the perturbation Lq. Thus the total 
Lagrangian is given by 


L = L c + Lq 


with components 


Lq = r d ~ l ( 2 - (Ed z E* - E*d z E) + |V £| 2 



and 

Lq = *r d_1 (W + e^ 4 ). 


( 8 ) 

(9) 


( 10 ) 


We want to point to the analysis presented in 1501153153 ] for comparable applications of the non-conservative La¬ 
grangian method with generalized forces in the context of Ginzburg-Landau type or complex Gross-Pitaevskii-type 
equations. We take <[8j) and apply Hamilton’s principle of vanishing variation, i.e. S(f /(Lq + L Q)dzdr) = 0, which 
is obtained for the minimizing condensate wave function (minimizer). The exact minimizer is a solution of the 
Euler-Lagrange equations, obtained from the Lagrangian L of the system by applying 



( 11 ) 


Although we do not know the exact form of the minimizer we can approximate the minimizer by an ansatz similar to 
the bright soliton solution for equilibrium condensates 0, the results will the later be compared to numerics yielding 
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the numerically generated minimizer justifying our ansatz. To proceed this way we first optimize over the transverse 
coordinates [5T1 1521 


* 


d ( d 


L- 


d 


drj \dE'* dE 


L = 2Re / dr 


dE 


8 


drj dE 1 * 


( 12 ) 


with i] = r, z where the r.h.s. corresponds to the non-conservative part EH • Furthermore by denoting L c = f dr L c 
we effectively consider EH 


d fdL, 


dz \ dv' 


dLc 

dv 


= 2 Re 


J r d 1 dr (iSE + ie\E\ 2 E) 


dE* 

dv 


(13) 


where the generalized forces on the r.h.s. represent the non-equilibrium processes perturbing the condensate. Here 
v = A(z),W(z),P(z),C(z) for the nonlinear bright soliton ansatz of the condensate wave function representing an 
approximate minimizer of (| 8 |) that is given by 


E(z, r) = A(z) ■ sech 


G V(z) 


exp (iP(z) + iC(z)r). 


(14) 


We insert this ansatz in L c and subsequently plug the expression for L c into the l.h.s of eq. |13| Then variation over 
A yields 


P' = 


-1 + 2A 2 W 2 - 3 C 2 W 2 
3W 2 


utilizing here the simple notation ' = d z . By variation regarding W we obtain 

1 


W = ± 


V-A 2 + 3 C 2 + 3P' 


(15) 


(16) 


and variation regarding P yields 


d z (A 2 W) = 2 A 2 SW+ -A 4 eW. 


(17) 


Finally variation regarding C gives us the formula 


4 A 2 CW = 0. 


(18) 


Expression (18) implies that C = 0 as long as we have a non vanishing amplitude and do not consider the homogeneous 
solution. We note that our results are consistent with [50:, [52] [53] and point out that the sister case of dark solitons 
in non-equilibrium BEC has been considered in [ 52] . For vanishing time dependence of the density profile, i.e. A and 
W are stationary, we obtain a simple formula for the amplitude of the bright soliton in terms of the linear pumping 
(reduced by the linear decay) 6 of polaritons/atoms into the condensate phase and the density dependent decay of 
strength e = — |e|, 


A = ± \W\- 


(19) 


One observes in (19) that the amplitude is higher for increased linear absolute gain S and lower for increased density 
dependent decay e. The simple limit <5 —^ 0, yields a vanishing occupation of the condensate phase, because A — l 0 
consistent with experiments. For the width of the bright soliton we obtain 


w = ±A = ± 

A 


4 N 
3 5 ’ 


( 20 ) 


which differs by a factor of \/2 from the corresponding dark soliton expression [52]. As the amplitude of the condensate 
increases, the width of the bright soliton decreases, due to this reciprocial relation. Finally the phase of the non¬ 
equilibrium bright soliton satisfies, 


A 2 

P{z) = -z- 


3 5 

const = - tZ - 

4 £ 


const. 


( 21 ) 


and increases quadratically with amplitude. 
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ATTRACTIVE NON-EQUILIBRIUM SPINOR CONDENSATE 

While we have introduced the non-conservative Lagrangian method for solitons in attractive BEC, next we consider 
the complete spinor system ([!]) for one spatial dimension [15] and formally proceed similarly as in 52]. Using the 
Keeling/Berloff non-equilibrium terms for the governing condensate i> we generalize them for the spinor condensate 
to 


('yn R - T) ~ i(5 + ea' 1 (|V’+| 2 + |V’-| 2 )), 


( 22 ) 


which is the mathematical form we shall analyst in the following. By applying analogous rescaling as in the coherent 
case ([ 6 ]) we obtain the spinor PDEs 


id z E ± = -AE ± - {\E ± \ 2 + 1 \E T \ 2 )E ± +i5E ± + ie(\E ± \ 2 + l2 \E T \ 2 )E ± , 


(23) 


fully describing state-of-the art polariton experiments in appropriate regimes and which are analogous to atomic BEC 
mean field models [2]. Here 7 = 02 / 0:1 is the magnitude of cross-scattering between the two spin components and 
72 models misalignment between self- and cross-spin decay. Correspondingly, the non-conservative Lagrangians (for 
dissipative and conservative part) for the attractive spinor condensate are 


Lc= E ( E ^ E *k - E ± d * E k) + |V£fc | 2 - ^ - l\E+\ 2 \E-\ 2 , 

fce{+,-} V J 

including cross-interaction terms in contrast to the coherent states formalism above, and 


L Q= E iUE k \ 2 +e^)+iei2\E + \ 2 \E_\ 

ke{+,-} V ' 


(24) 


(25) 


Here the last term in (251 is due to interspecific competition between the two spin components. We adapt ansatz (14) 
to the spinor case, writing 4/' = (E + ,E-) we define its two components by 


E±(z,t) = A±(z) • exp - 


{r ± r 0 f 
W ± (z ) 


exp(zP ± (z) + iC±{z)r), 


(26) 


where we model the bright solitons in terms of a gaussian instead of the sech(a;) discussed in the previous section. 
The half-bright solitons are parametrized through the time dependent amplitude A ±, width W ±, phase P ±, spatially 
dependent phase factor C± and a distance between the half-solitons given by r$. This ansatz describes two half-bright 
solitons positioned at r = ±r 0 correspondingly. To integrate the cross-interaction term we proceed as follows. First 
we obtain the integrated conservative part of the Lagrangian as in the previous section. The result is 


Lc= E ^ A lJ^r(- A lW, + 2V2(l + W a (C 2 + P^))-' 1 AlA 2 _ 



Sr- 2 


W a + W w 


(27) 


CT — + ,— 

When the distance between the bright solitons is tq —> 00 the cross-interaction is effectively zero, which is plausible 
due to finite sized widths of the bright solitons. Similarly as 7 —> 0 (27) uncouples as well the conservative part of 
the two components and we consider the case of two independent half-solitons, which resemble the single component 
case. In accordance with the treatment of the previous section the generalized integrated Lagrangian equations are 


d fdL c 
dz I dv' 


dLc 

dv 


= 2Re / r dr[iSE^ + ie(\E a \ + 7 2 |£W| )E 


dE* 

dv 


(28) 


We remark that to obtain complete decoupling between both spin components we consider the limit 7 —► 0 and 
72 —> 0. Here 72 —» 0 corresponds to the case of vanishing interspecific competition 152] , In general (for 7 > 0 and 
72 > 0 ) component + will impact the parameters of the bright soliton within the other component — and vice versa, 
i.e. the system is coupled through cross-interaction and cross-competition. 

Using 


(28) and the 8 parameter ansatz (26) we obtain a set of 8 coupled ordinary differential equations for 
those parameters: Variation over W± yields 
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. Al 

P' — _±_ 

* 2y/2 


-Cl 


1 

w± 




w± 


16 exp ^ w± + ° Wt ) to \Jw± + w T ^ i 7 " 


^(Wi + W^ 


(W± 


W T ) 3 


(29) 


with 7 denoting the cross-interaction strength, and variation of A± gives 


C± = ±\ 


V2A 2 ± -2Pl - AC + 7exp (- 1 jg^ )- 


9 42 / 1 

I/ W± 


VU ± ' 


Finally variation of P± and C± yields correspondingly 

d z 


(AlVW)=Al ( Sy /wZ + ‘U y /wZ + J*JML 
V v V v w± + w. 


: exp - 


8 r 0 2 


where 72 denotes the interspecific competition strength, and 

roV2nW± 


A 2 ± C ± V2^VWl = E12A 2 ± A 2 - 


^—^ {W± + w T ) eXP \ W± + W T 


w± + w* 


8 r 2 0 



(30) 


(31) 


(32) 


Let us now derive a stationary density solution, where d z A± = 0 = d z W±. Assuming ro = 0 the r.h.s of (32) vanishes, 
so C± = 0 for a nontrivial non-homogeneous solution, which reduces the system to 6 free parameters. Furthermore 
we suppose W± = W T . Using (311 we obtain for the amplitudes 


^± = 


y/26 ( -1 


I + 72 


(33) 


Note that e = —\e\. Formula (331 implies that the higher the interspecific cross-competition, the lower the amplitude 
of each component , while the dependence on 6 and e corresponds with the findings for the spin coherent case. The 
widths of the two bright solitons satisfy 


W± = 


4e /1 + 72 
<5 V 1 + 7 


(34) 


which again is in accordance with the findings made in the previous section. However here for larger inter and intra 
component decay the solitons widen and as pumping increases the bright solitons narrow. Finally the phases are given 

by 


A 2 (l + 7 ) S 

P±{z) =-—- + c 0 = - — 


4\p2 


7+1 


4e V1 + 72 


Cq, 


(35) 


where Cq is the constant of integration. Note that inserting these expressions into (26 1 yields the two condensate wave 
functions of the spinor non-equilibrium system. 

In the following section we show numerical phenomenology of bright solitons and spinor bright solitons within 
non-equilibrium condensates and compare some numerics with the above analytical results. 


NUMERICAL RESULTS 


In this section we present direct numerical solution to Eqs.(2[3) that were computed with an adaptative step 
Runge-Kutta-based solver. We worked here with a spatial resolution of h = 0.25^m. Let us first recall the case of a 
conservative closed condensate (e.g. an atomic condensate without exchange of particles or a polariton condensate 
with very long lifetimes due to high quality cavities). This setting corresponds to the assumptions P = 0 and T = 0 
in Eqs. In the coherent component case, which requires to set in addition 012 = 0, the systems reduces to 
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FIG. 1: Bright solitons in an equilibrium scalar condensate. The colormap shows the chemical potential g{x,t) = Qi n(x,t) 
(meV). (a) Stable soliton solution from i/V (b) Unstable breathing obtained imposing a too large intial condition, (c) Bright 
soliton attraction and oscillations. 


considering the single component evolution, since 'i/> + = ?/>_ = il>. It is well-known that bright solitons or rogue waves 
naturally emerge in attractive condensates from the collapse of an inhomogeneous particle distribution |9j. We follow 
this study and let us consider in the remainder of this work typical exciton-polariton condensates parameters to be as 
explicit as possible and without loss of the generality for the observations made here. In these systems the temporal, 
spatial and energy scales are respectively a few ps, fim and meV s. Typical scenarios of different development in 
time of the bright soliton is shown in Fig. |TJ panel (a) shows the unchanged form of the wave function evolved as 
described above starting from an initial bright soliton wavefunction ?/{, (x) = -y/riosech(x/£) where £ = h/y/rnfi is the 
healing length setting the soliton size and p = aqno is the chemical potential. Panel (b) shows a realization where 
the initial wavefunction is set twice the size, i.e. with a spatial extension of 2£ ~ 2/un. Here the interactions are 
not properly compensated by dispersion and the wavepacket/bright soliton starts collapsing. Subsequently it expands 
again when dispersion becomes dominant forming a clear breathing pattern that demonstrates their intrinsic instability 
due to the attractive interactions. Around 120 ps after starting the process the finite numerical tolerance is enough 
to spontaneously break the symmetry. Then the bright soliton escapes to w.l.o.g. the left transforming the excess 
of energy into motion. Finally in the panel (c) we have initially prepared two bright solitons separated by 4 £ and 
we observe their expected attraction. Indeed density maxima are bound to attract each other when a± < 0. The 
attraction infers a velocity to the solitons that cross each other and repel again performing dipolar oscillations until 
they split at t ~ 100 ps. 

Now we consider the equilibrium spinor case where both components are taken into account. In correspondence 
with above findings, in the case of negligible cross interactions, where = 0 a stable spatial soliton in both component 
can be set via ip± = i/b as shown in Fig. [2] (a). Here the left (right) part shows the cr + (<r_) component dynamics. 
Now if the solitons in each component (half-solitons) are spatially separated, they attract for u-i < 0 and repel when 
0.2 > 0 as shown in Fig[2] (b) and (c) respectively as opposed to the dark soliton case m- First note that in contrast 
to the dark soliton case where the effective magnetic held was shown to accelerate the half-solitons appearing on top 
of a nonzero density background, here it simply results in a pseudospin precession namely a weak population exchange 
between the two components. Secondly note that the systems chemical potential is modified as p = Uo(ai +a2)/2 in 
the spinor case. 

These numerics provide us intuition for the non-equilibrium case and we can already imagine that the competition 
between the intra-component repulsion (Fig. [lj and cross interactions (Fig. [2]) combined with the breathing effect 















FIG. 2: Bright solitons in an equilibrium spinor condensate, the colormap shows the chemical potential, (a) Stable soliton 
solution from = ipb- Half-solitons (a) attraction for ai < 0 and (c) repulsion for <22 > 0. 


can lead to rich dynamics of the non-equilibrium spin sensitive BE condensate. 

Let us now treat the full problem as described by Eq. [2] including the effects of the particles being scattered into the 
condensate modes, i.e. formally P / 0 and decay r ^ 0. The initial population profile is set by the pump amplitude 
fixed to A P = 2501^. A stable soliton solution can be obtained from a narrow pump spot where crp = 2 /im as shown 
in Fig. [ 3 ] (a) where overlapping stationary bright solitons are formed after 100 ps in both components. Interestingly 
the sustaining pump can wash out small perturbations in such a way that the soliton stability region is narrow but 
not singular. The inset shows a slice at t = 500ps (white line) with the analytical sech (dashed-red line) from which 
we deduce a soliton width of W = 0.85 jirri. 

Supposing the reservoir favors one component over the other, as it would be the case e.g. when pumping the 
condensate with a spin polarized pump laser or spontaneously in trapped condensates far away from the pumping 
spots 1221, one can introduce a symmetry breaking between the two spin components. In Fig. [ 3 ] (b) we have increased 
the transfer rates 7 to the a + component by 1% and considered a 2 > 0. This scenario results in an imbalance between 
the a± polariton/atoms populations and not only the er + component does not match the stable bright soliton anymore 
leading to dominant dispersion but also it is reduced by the cross repulsion brought by <22 >0. As a result the <r+ 
population vanishes leaving a stable half-soliton in the er_ component of width W = 0.55 /im (see inset), which is in 
accordance with the experimental findings in m ■ However in l27| spin selection happened to be spontaneous without 
any imbalance of the pump rates as discussed here. In the case <22 < 0 (not shown) the cross attraction prevents the 
spreading to occur leading to the overlapping stationary solitons. 

Finally in the most common case where the pump spot size does not match the the optimal soliton condition as 
illustrated in Fig. [ 3 ] (c) and (d), we obtain a very complex fireworks dynamics resulting from the multiple cross/self¬ 
interactions and competition and the breathing effects. Here we used the parameters up = 15 [im for the spot size in 
the case <22 < 0 in panel (c) and <22 > 0 in panel (d). We note that including a weak noise to the system dynamics, 
that is present in state-of-the-art experiments, leads to random patterns of the bright solitons, which can be utilized 
for the realization of a true random number generator. 


CONCLUSIONS 

We have presented an approximate analytical treatment for the parameters of spin sensitive non-equilibrium bright 
solitons that is based on a non-conservative Lagrangian formalism. This Lagrangian formalism has been extended 
from a previous work to include the attractive interaction case of the non-equilibrium spinor condensate field and 
turned out to provide reliable results. Explicit ordinary parameter differential equations for the amplitudes, widths 
and phases were derived for the coherent and spin sensitive case and interpreted in terms of key system parameters. 
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FIG. 3: Driven dissipative solitons, the colormap shows the chemical potential. The a± label refers to the ± component, (a) 
Stable soliton solution obtained with ap = (b) Stable half-soliton for 0:2 < 0. Unstable dynamics for crp = 15/rm. (c) 

a 2 < 0, (d) 02 > 0. The insets in panel (a) and (b) shows a slice at t = 500 ps (white line) fitted by the sech analytical function 
(dashed-red line). In both cases the mean absolute error between the numerically and analytically generated curves is less than 
10~ 5 meV. 


Particularly explicit expressions for a stationary solution of the bright soliton in a non-equilibrium environment were 
stated. In addition the scenario of non-equilibrium half-bright solitons was elucidated analytically and we found 
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corresponding parameter equations and explicit static spinor expressions that include the coherent case as limit for 
vanishing cross-interactions and competition. The explicit solutions allowed us to clarify the role of cross-interactions 
and interspecific competition between the two spin components of the condensate for the system parameters, the 
latter inducing a decrease of amplitudes. Furthermore we have demonstrated numerically the generation and high 
control over bright soliton emergence in non-equilibrium systems and highlighted several effects important for the 
fate of the bright solitons. The existence of stationary equilibrium bright solitons in the spin coherent case as well as 
breathing solutions with subsequent spontaneous symmetry breaking was shown, which is a simple manifestation of a 
true random number generator. In addition we considered the non-equilibrium dynamics of spinor condensates, where 
we have found for instance the novel phenomenon of bright soliton fireworks dynamics. Finally we have compared 
the analytical and numerical expressions and obtained a very good match of the presented results. 
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